numerical solution of incompressible boussinesq equations using fourth-order compact scheme: lock exchange flow

Authors

سرمد قادر

استادیار، گروه فیزیک فضا، مؤسسة ژئوفیزیک، دانشگاه تهران، ایران ابوذر قاسمی ورنامخواستی

دانش آموخته کارشناسی ارشد، گروه فیزیک دریا، دانشکده علوم دریایی، دانشگاه تربیت مدرس و محقق پژوهشگاه هواشناسی و علوم جو، نور، ایران محمدرضا بنازاده ماهانی

استادیار گروه فیزیک دریا، دانشکده علوم دریایی دانشگاه تربیت مدرس، نور، ایران داریوش منصوری

مربی گروه فیزیک دریا، دانشکده علوم دریایی دانشگاه تربیت مدرس، نور، ایران

abstract

in recent years, the number of research works devoted to applying the highly accurate numerical schemes, in particular compact finite difference schemes, to numerical simulation of complex flow fields with multi-scale structures, is increasing. the use of compact finite-difference schemes are the simple and powerful ways to reach the objectives of high accuracy and low computational cost. compact schemes, compared with the traditional explicit finite difference schemes of the same order, have proved to be significantly more accurate along with the benefit of using smaller stencil sizes, which can be essential in treating non-periodic boundary conditions. applications of some families of the compact schemes to spatial differencing of some idealized models of the atmosphere and oceans, show that the compact finite difference schemes are promising methods for numerical simulation of the atmosphere–ocean dynamics. this work is devoted to the application of a fourth-order compact finite difference scheme to numerical solution of gravity current. the governing equations used to perform the numerical simulation are the two dimensional incompressible boussinesq equations. the two-dimensional lock-exchange flow configuration is used to conduct the numerical simulation of the boussinesq equations. the lock-exchange flow is a prototype problem which has been studied numerically and experimentally by many researchers. for the spatial differencing of the governing equations the second-order central and the fourth-order compact finite difference schemes are used. the predictor-corrector leapfrog scheme is used to advance the boussinesq equations in time. the boundary condition formulation required to generate stable numerical solutions without degrading the global accuracy of the computations, is also presented. the fourth-order compact scheme is compared in detail with the conventional second-order central finite difference method. qualitative comparison of the results of the present work with published results for the planar lock-exchange flow indicates the validity and accuracy of the fourth-order compact scheme for numerical solution of the two-dimensional incompressible boussinesq equations.

Upgrade to premium to download articles

Sign up to access the full text

Already have an account?login

similar resources

A Fourth Order Scheme for Incompressible Boussinesq Equations

1 Institute for Physical Science and Technology and Department of Mathematics, University of Maryland, College Park, Maryland 20742. 2 Institute for Scientific Computing and Applied Mathematics and Department of Mathematics, Indiana University, Bloomington, Indiana 47405. E-mail: [email protected] 3 Department of Mathematics, University of Michigan, Ann Arbor, Michigan 48109. Received February ...

full text

numerical solution of unsteady and non-linear rossby adjustment problem using fourth-order compact maccormack scheme

the compact finite difference schemes have been found to give simple ways of reaching the objectives of high accuracy and low computational cost. during the past two decades, the compact schemes have been used extensively for numerical simulation of various fluid dynamics problems. these methods have also been applied for numerical solution of some prototype geophysical fluid dynamics problems ...

full text

A Compact Fourth - Order Finite Difference Scheme for the Steady Incompressible Navier - Stokes Equations

We note in this study that the Navier-Stokes equations, when expressed in streamfunction-vorticity fonn, can be approximated to fourth--order accuracy with stencils extending only over a 3 x 3 square of points. The key advantage of the new compact fourth-order scheme is that it allows direct iteration for low~to-mediwn Reynolds numbers. Numerical solutions are obtained for the model problem of ...

full text

A High Order Numerical Scheme for Incompressible Navier-Stokes Equations

To solve the incompressible Navier-Stokes equations in a generalized coordinate system, a high order solver is presented. An exact projection method/fractional-step scheme is used in this study. Convective terms of the Navier-Stokes (N-S) equations are solved using fifthorder WENO spatial operators, and for the diffusion terms, a sixthorder compact central difference scheme is employed. The thi...

full text

Analysis of a fourth order finite difference method for the incompressible Boussinesq equations

Cheng Wang, Jian-Guo Liu, Hans Johnston 1 Department of Mathematics, University of Tennessee, Knoxville, TN 37996, USA; e-mail: [email protected] 2 Institute for Physical Science and Technology and Department of Mathematics, University of Maryland, College Park, MD 20742, USA; e-mail: [email protected] 3 Department of Mathematics, University of Michigan, Ann Arbor, MI 48109, USA; e-mail: hansjo...

full text

My Resources

Save resource for easier access later


Journal title:
فیزیک زمین و فضا

جلد ۳۷، شماره ۱، صفحات ۰-۰

Keywords
in recent years the number of research works devoted to applying the highly accurate numerical schemes in particular compact finite difference schemes to numerical simulation of complex flow fields with multi scale structures is increasing. the use of compact finite difference schemes are the simple and powerful ways to reach the objectives of high accuracy and low computational cost. compact schemes compared with the traditional explicit finite difference schemes of the same order have proved to be significantly more accurate along with the benefit of using smaller stencil sizes which can be essential in treating non periodic boundary conditions. applications of some families of the compact schemes to spatial differencing of some idealized models of the atmosphere and oceans show that the compact finite difference schemes are promising methods for numerical simulation of the atmosphere–ocean dynamics. this work is devoted to the application of a fourth order compact finite difference scheme to numerical solution of gravity current. the governing equations used to perform the numerical simulation are the two dimensional incompressible boussinesq equations. the two dimensional lock exchange flow configuration is used to conduct the numerical simulation of the boussinesq equations. the lock exchange flow is a prototype problem which has been studied numerically and experimentally by many researchers. for the spatial differencing of the governing equations the second order central and the fourth order compact finite difference schemes are used. the predictor corrector leapfrog scheme is used to advance the boussinesq equations in time. the boundary condition formulation required to generate stable numerical solutions without degrading the global accuracy of the computations is also presented. the fourth order compact scheme is compared in detail with the conventional second order central finite difference method. qualitative comparison of the results of the present work with published results for the planar lock exchange flow indicates the validity and accuracy of the fourth order compact scheme for numerical solution of the two dimensional incompressible boussinesq equations.

Hosted on Doprax cloud platform doprax.com

copyright © 2015-2023